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Abstract 

We investigate the decoupling limit in the DGP model of gravity by study- 
ing its nonlinear equations of motion. We show that, unlike 4D massive grav- 
ity, the limiting theory does not reduce to a sigma model of a single scalar field: 
Non-linear mixing terms of the scalar with a tensor also survive. Because of 
these terms physics of DGP is different from that of the scalar sigma model. 
We show that the static spherically-symmetric solution of the scalar model 
found in|hep-th/0404159 , is not a solution of the full set of nonlinear equations. 
As a consequence of this, the interesting result on hidden superluminality un- 
covered recently in the scalar model in hep-th/0602178, is not applicable to 
the DGP model of gravity. While the sigma model violates positivity con- 
straints imposed by analyticity and the Froissart bound, the latter cannot 
be applied here because of the long-range tensor interactions that survive in 
the decoupling limit. We discuss further the properties of the Schwarzschild 
solution that exhibits the gravitational mass-screening phenomenon. 



1 Introduction and summary 



In a simplest perturbative expansion non-linear interactions in the DGP model PP 
become strong at a certain scale determined by the graviton mass, coupling constant 
and physical parameters of a problem at hand (e.g., mass/energy of the source) |2]. 
The strong interactions and its scale can be seen in a breakdown of the classical 
perturbation theory around a static source (in analogy with a similar phenomenon 
observed in Ref. in the context of massive Fierz-Pauli gravity (FP) |3|). More 
generically, this can be understood in terms of non-linear Feynman diagrams that are 
enhanced by inverse powers on the graviton mass [2^ . This perturbative breakdown, 
and classical resummation of the corresponding large diagrams, is what makes the 
model agreable with predictions of General Relativity j2], 0-0 , and yields tiny but 
potentially measurable deviations from it [IHl IHl 112] , [inilll]- 

It is instructive to consider the decoupling limit of the DGP model, as it was 
done in Refs. |T3] and [1^] , motivated by the original studies of the decoupling limit 
of massive gravity in Ref. |17j. It was argued that in this limit the DGP model 
reduces to a non-linear scalar model, while the self-interactions of a tensor mode 
vanish fT^ ITT)] . The remaining equation for the scalar field reads: 
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where A = Mpim^ is held fixed while Mpi —>■ oo, rric 0, and rric plays the role of the 
graviton mass (the above equation is given here in a source- free region). Although, 
the results of ^3 E] are similar in spirit to those found in in the context of 
massive gravity, there are also important differences between the scalar models of 
Refs. |T7] and 'T31C5]- These differences were shown to reproduce [THIITU] the known 
non-linear ghost-like instability in the FP theory [20], |2I]; and were argued to be 
responsible for its absence in ((T]) [rS] HH] • 

In spite of all the above, the formal similarity between the two theories seems 
somewhat surprising because of fundamental differences in their Lagrangians (see 
more on this in section 6). Most importantly, however, there exist an expression for 
the metric J2] such that: (a) It is an exact solution of the nonlinear DGP equations 
on and near the brane; (b) It is a perfectly regular solution in the decoupling limit 
(see section 5); (c) Yet, it does not satisfy Eq. ([T]). 

The above observation motivated us to re-examine the decoupling limit of the 
DGP model in the present work. By analyzing the full set of nonlinear equations 
we will show below that one of the equations that survives in the decoupling limit 
is more general than (P), and looks as follows: 
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Here R = MpiR and R^i, = MpiR^^, where R and R^i, are the Ricci scalar and 
tensor of a spin-2 field. Eq. (j2J) reduces to Eq. ((TJ) if, e.g., both Ricci tensor and 
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scalar vanish. However, this does not have to be the will be discussed below. 

For instance, the Schwarzschild solution of [121 Ej gives R and components of R^y 
that are nonzero even in the decoupling limit, and satisfies (0), but not ((T)). 

The issue of small fluctuations about classical backgrounds is different in (|T)) 
and (j21). Eq. ((T)) implies that the fluctuations of the vr field are decoupled from 
those of the tensor and carry an independent physical meaning. On the other hand, 
Eq. Q shows that the fluctuations of the vr field mix with those of a tensor. In 
fact, as we will see below, the vr field can be completely absorbed into a tensor 
field, which survives as a nonlinearly interacting mode and propagates five physical 
polarizations. 

More importantly, however, there exist yet another non-linear equation that 
survives in the decoupling limit. The latter is somewhat involved to be presented 
here without lengthy explanations (see Eq. (fTH|) ). Consistent solutions should satisfy 
that equation too. We will show that the spherically-symmetric solution of Eq. (0) 
found in Ref. |^, although obeys does not satisfy Eq. (fTH|l . Hence, it is not a 
solution of the DGP model^. On the other hand, the solution of TH" satisfies both 
© and (HHD- 

As was recently shown in Ref. [221 ? the scalar model (P) presents a very interest- 
ing field-theoretic example: In spite of its appearance as a local, Lorentz-invariant 
effective field theory, it exhibits certain hidden non- locality [221 • The latter can 
manifest itself in superluminal propagation of perturbations on the background de- 
fined by the spherically-symmetric solution of [in]. If the DGP model were to reduce 
in the decoupling limit to the scalar theory (0), then, any experimental test of this 
model would present a window of opportunity for discovering small effects of an 
0{1) superluminality in gravitation [22] ■ 

While this intriguing possibility can exist in the scalar sigma model [22] 5 our 
results show that it is not offered by DGP. The solution of [TEj, and its fluctuations 
that manifest superluminal propagation, are not consistent solutions of the DGP 
model. More generally, the scalar model violates positivity constraints [22] imposed 
by analyticity and the Froissart bound, however, the latter cannot be applied to the 
present model of gravity because of the long-range tensor interactions that survive 
in the decoupling limit. 

The paper is organized as follows. In section 2 we discuss the ADM formalism 
for DGP and write down explicitly some of the key nonlinear equations. In section 3 
we take the decoupling limit of those equations and derive Eq. In section 4 we 
deal with the remaining nonlinear equations of the model and find their decoupling 
hmit. We also show that the solution of Ref. [T^] does not satisfy these equations. 
In section 5 we discuss the physical interpretation of the Schwarzschild solution of 
Refs. [121 m] in the decoupling limit. We emphasize that this solution, although 

^The solution of [Iftj . in our nomenclature, is a solution of the equations of motion that are 
obtained from the sigma model Lagrangian in which the tt field is decoupled from the tensor field 
(i.e., there are no mixings between them) |15[I16| . In the solution discussed in W the tt and tensor 
fields are sourcing each other. 
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naively seems counterintuitive, actually has a clear physical interpretation in terms 
of the screening effects. In section 6, we discuss some interesting open questions. 



2 DGP in the ADM formalism 

Let us start with the action of the DGP model |I] in the ADM formalism (see, 

e.g.. El, my- 

S=^Jd'x dy^g [r 6{y) + '^n{r + K'- K.^K^'')^ . (3) 

Here, the (4 + 1) coordinates are x^^ = (x^, y), /i = 0, . . . , 3; g and R are the 
determinant and 4D curvature for the 4D components g^j,y{x,y) of the 5D metric 
gAB{x,y). K = g^^Kfj_y is the trace of the extrinsic curvature tensor defined as 
follows: 

K,, = ^ {dyg^u - V^iV, - V.iV^) , (4) 

and is a covariant derivative w.r.t. the metric g^^j. We introduced the lapse 
scalar field A^, and the shift vector field A*"^ according to the standard rules: 

g^y = N, = g.^N" , gyy = + g,, ^ N" . (5) 

Equations of motion of the theory are obtained by varying the action w.r.t. g^j^,^, 
N and A^. Here we start with a subset of two equations, the junction condition 
across the brane, and the {yy} equation that can be obtained by varying the action 
w.r.t. A^. The former reads as follows: 

Gf,u - mc{K^^ - g^uK) = T^^/Mp^ , (6) 

where is the 4D Einstein tensor of the induced metric gfj,^ and T^j, is the matter 
stress tensor. The {yy} equation takes the form 

R=K^-Kl. (7) 

Note that © is valid only at y = while ([7j) should be fulfilled for arbitrary y. The 
terms with the extrinsic curvature contain derivatives w.r.t. the extra coordinate 
as well as the yy and fiy components of the metric. Nevertheless, one can deduce a 
single equation that contains the 4D induced metric only! This is done by expressing 
from (jH)) K and Kf^i, in terms of R and R^^i, outside of the source 

K,^ = ^^^^^^^, K=^, (8) 
rrir orrir 



and substituting these expressions into ((Zj). The result at y = 0^ is: 

R^ - 3R' 

We will study the decoupling limit of this equation below 
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3 (De) coupling limit 



The purpose of this section is to show that the nonhnear interactions of the tensor 
field in the DGP model do not disappear in the decoupling limit. 
The decoupling limit in DGP is defined as follows |T3j: 

Mpi-^oo, rUc^O, M^cx), (10) 

where M is a mass of a source entering the stress-tensor T and the following quan- 
tities are held finite and fixed 

A^{M,,miy/\ (11) 
To take this limit in Q, we multiply both sides of that equation by Mpi and obtain 

R = (12) 

where we defined R = Mp\R and R^i, = M-pxR^^. Now we are ready to proceed to 
fllOpilj) . Before doing so, it is instructive to make parallels with FP gravity. The 
analog of the mass term in the DGP action (jH)) is the one containing squares of the 
extrinsic curvature. Looking at the expression for the extrinsic curvature (0} one can 
think of g^y in DGP to be an analog of the tensor field of FP gravity in which the 
Stiickelberg field is manifestly exposed in the mass term ^7] . Then, is an analog 
of the vector-like Stiickelberg field as far as the reparametrizations with the gauge- 
function C^(x,y) = (C^(a;, y), 0) are concerned. The longitudinal component of A^"^ 
should play the role similar to that of the longitudinal component of the vector-like 
Stiickelberg field of FP gravity (note that there is an additional Stiickelberg field N 
in DGP, the important role of which will be discussed in section 6). 

Equation (HH), on the other hand, does not contain any of the Stiickelbeg fields 
but the tensor g^^ which in the decoupling limit has a conventional scaling^. That 
is why it is straightforward to take the limit directly in (|12p. For this we recall that 

(hUh\ 

Rf" = \ - B^d'K. - duSfKf + d^dS) + O , (13) 

where we defined a field /i^j^ = M-p\{g^j_y—ri^y), which has the canonical dimensionality 
and is held fixed in the limit^. 



^If we were to do perturbative calculations, this would imply the gauge choice for which the 
Stiickelberg fields are manifestly present. 

■^The signs in the expressions (|13ll and the definition of given above determine our convention 
for the sign of the curvature tensor. 
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Note that we have not done any small field approximation, but merely expanded 
the nonlinear expressions for R and Rfj_y in powers of h and took the decoupling 
limit. Because all the nonlinear terms in R and R^y are suppressed by extra powers 
of Mpi, the curvatures reduce to their linearized form. 

The expressions in flT^ . should be substituted into (fT^ . To make closer contact 
with the TT language of Section 1 we perform the following shift of the h^j^ field 

h^,y = h^u + V^iv-T^ ^ (14) 

where the tt field has a canonical dimensionality and is held fixed in the decoupling 
limit. With this substitution Eq. (fT^ turns into ^ which we repeat here for 
convenience 

3°- + ^= + 3A3 + • ^^^^ 

R and R^u denote the linear terms on the right-hand-sides of (fT3|) where h^y is 
replaced by h^y. 

Eq. (|T5|) shows that a tensor field would have non-linear interactions in the 
decoupling limit as long as there are no additional equations constraining all the 
components of R^^y to be zero. We will turn to the remaining equations in the next 
section and show that they do not necessarily imply vanishing R^y. Before that we 
would like to make a comment on small fiuctuations about classical solutions. The 
TT field in (fT3j) can be reabsorbed back into the tensor field h^y using ((T3j). Small 
perturbations on any background in the decoupling limit are those of hp,y. The given 
background is what defines the light-cone, and there are no additional degrees of 
freedom that could propagate outside of that light-cone. 



4 More bulk equations 

Some of the bulk equations have not been considered in our discussions so far. 
These are the {/iy} and bulk {/iz/} equations (the {^iv} equation on the brane (jH} 
has already been taken into account). We will discuss them in the present section. 
We start with the {^y} equation which for arbitrary y reads as follows: 

V,K = V'K^y. (16) 

The covariant derivative in the above equation is the one for g^y. At ?/ = Eq. (fTBj) 
is trivially satisfied due to (jH)). For y 0, (fTBj) sets the relation between N^, N and 
Qf^y. Hence, (fT^ gives a relation between these quantities for both y = and y 0. 

One can use the bulk {yy} equation ((Tj) to determine N in terms of A*"^ and 
gfj_y, and then use (fTHj) in order to express in terms of g^y. If so, there must 
exist one more equation which should allow to determine the bulk g^j^y itself. This 
is the bulk {fiiy} equation, to which we turn now. The latter can be written in a 
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few different ways. For the case at hand, the formahsm by Shiromizu, Maeda and 
Sasaki |2n| is most convenient. In this approach, the bulk {/xi^} equation and the 
junction condition can be combined to yield a {/xz^} equation at y ^ 0~^. This gives 
a "projection" of the bulk {yuz/} equation onto the brane. Since the bulk itself is 
empty in our case, this is the most restrictive form one can work with. The equation 
reads as follows |2B|: 

= KK^y - KP^Ky, - \g,.{K^ - K^,) - E,, , (17) 

where all the quantities are taken at ?/ = 0^. G^^, denotes the 4D Einstein tensor of 
the metric g^^, E^^ denotes the electric part of the bulk Weyl tensor, projected onto 
the brane. An important property of the Weyl tensor is that it is invariant under 
the conformal transformations. Moreover, Efj,^ is traceless. 

We now turn to the decoupling limit in this equation. For this we multiply both 
sides of (fTTj) by Mpi, exclude K and K^^ using (jS)), and take the limit (fTUl ITT|) 
holding the canonically normalized components of g^^ fixed. The resulting equation 
reads: 

Gf,u = - Ef,„ , (18) 

where 

1-- 1- - 1 -2 

Bfiu = —G^aG'^ + -G'^Gfj^i, + -Tj^^GafsG'^ — -Tlf_iu{G'^) , (19) 

and Gf^i, = M-piG^y, E^^ = M-p\E^y. This is the equation that has to be satisfied by 
any consistent solution. 

It is straightforward to check that the trace of the above equation is nothing but 
fll2|l . Moreover, as in Eq. (fT^ . there are nonlinear interactions of the tensor field 
that survive in the decoupling limit in (jlSp . This is true irrespective of the explicit 
form of E, as long as it's traceless. Once the values of and Ny are determined, 
as described at the beginning of this section, one can calculate E^^ in terms of g^u- 
Thus, Eq. (jl8|) turns into a single equation for the determination of g^^,. 

The above considerations show that the DGP model in general does not reduce 
to the scalar theory ((H) in the limit pOlll|) . This can also be deduced by looking at 
the solutions of the scalar model (Q) which are not solutions of ()18|). One example of 
this is the spherically-symmetric static solution of ((T)) found in jTH]. We will discuss 
this solution in the remaining part of this section. The ansatz corresponding to the 
solution of (TO] reads: 

^ _ _ 

Rt,u{h) = 0, G^^,{h) = d^duTi - T]f,^nn, K^^ = -^d^d^n . (20) 

One can calculate the expression for on the ansatz {E^i, comes out to be 
non-zero), and use that expression in (fTH|) . The resulting equation, after substitution 
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(H}, reduces to Eq. presented below. These calculations are tedious, and instead 
of them, we present here relatively easier calculations that produce the same result 
(this also serves as a self-consistency check of our result). 
Since the solution (j^UI) is written in terms of the extrinsic curvature, it is easier 
to check its compatibility with the bulk {/it^} equation also written in terms of this 
quantity. This equation takes the form 



G^u = \g,. {k' - Kip) + 2 (K'^K,^ - KK, 



V^V,N-g^,V^N dyi^iK^" -Kg^")) 



N + ^^"^'^ AV^^ • ^^^^ 

To check whether the solution of fT^ is consistent with ()21|) . we take the decoupling 
limit in the above equation with the substitution ^3 Il6j 

Ky = -'^d,Ii, hyy = -'^^dyli, H ( X , ^ ) = 6 " ^ ^^7^ ( X ) . (22) 

This can be expressed in terms of the lapse scalar and shift vector in the decoupling 
limit 

2 

Tfl 

N, = h,y, = gyy + ^{d,Ur . (23) 

The decoupling limit for gyy has been specified in Ref. ^H] only in the leading 
approximation gyy = 1 — 2m,cdyll / . However, for consistency, the subleading term 
in (PT|) is also needed, otherwise the ansatz itself would be inconsistent as it would 
violate some sacred properties of the bulk equation (j2H) (for instance, a non-linearly 
incomplete ansatz, after substitution into the bulk equation, does not respect the 
Bianchi identities and would not reproduce correctly the trace equation). On the 
other hand, in the next-to-leading order there is a unique non-linear completion for 
gyy that is consistent with the Bianchi identities and trace equation. The latter 
reads gyy = 1 — 2mcdyIl/A^ + ml{dyliy / . For this ansatz, the expression for 
in the decoupling limit reads as follows: A^ ~ 1 — nicdyll/A^ + m^((9^n)^/2A^. The 
latter combined with the other components above gives a consistent ansatz on which 
Eq. (PT|) looks as follows: 

d^,d,7i - rj^^OTi = -r]^, + ^ . (24) 

As before, trace of ()24|) gives precisely (0), i.e., the former equation is more general 
than (^. Hence, not all the solutions of satisfy (j23)- For instance, consider the 

■*Note that our curvature sign conventions are different from those of Ref. pS]. 
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{00} component of Eq. ()24|) and look at a static spherically-symmetric solution. 
Eq. reduces to 

2A3 ' ^ ^ 

where A denotes a 3D Laplacian and i,j = 1,2,3. The above equations is incom- 
patible with the time independent part of (P), unless both left and right hand sides 
of and are zero. The solution found in [16j does not satisfy this constraint, 
and, therefore, it does not satisfy the {00} component of Eq. (j^ . 

Likewise, one can show that the solution of ^H] does not satisfy the {ij} com- 
ponents of (|24p either. To demonstrate this, we consider the /i 7^ z/ part of Eq. 
(j24p and look at a spherically-symmetric static solution. Using the notations of ^H] , 
Ej = djTT = rjE/r, where j = 1, 2, 3, equation with i ^ j gives 

= . (26) 

The only nontrivial solution of ()2fi|l reads: E = cr, where c can be an arbitrary 
constant (including c = A^), which is not the solution of ^Hj- This is because the 
ansatz ()22j) is a pure gauge in the bulk only in the linearized approximation, but 
not in a full non-linear theory. In the next section we discuss a solution that exactly 
satisfies the above equations on and near the brane [131 El • 




5 Schwarzschild solution 

The Schwarzschild solution found in [T3] satisfies all the above nonlinear equations 
on and near the brane. Here, we will discuss some of the main features of this 
solution in the decoupling limit and give its further physical interpretation. 

For a static point-like source. Too = —M6^{x) with Tij = 0. A new physical 
scale emerges in this problem as a combination of Vc and tm 12] {i^m = 2GnM is the 
Schwarzschild radius and the Newton constant) 

n^iTMriy/'. (27) 

(This is similar to the Vainshtein scale in massive gravity P|). It is straightforward 
to check that this scale is finite and fixed in the decoupling limit (jlOlll|) . 

The metric on/near the brane was found exactly in Refs. P2til4j. Here, for our 
purposes it suffices to concentrate on the Newton potential alone 0(r) = /ioo/2. In 
the notations adopted in the previous sections 

goo = -l + hoo = -l + hoo/Mpi. (28) 

The exact expression for ^00 [IHI is a solution of Q- Furthermore, (|T7|) can be used 
to determine the off-diagonal and {yy} components of the metric near the brane. 
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as was done in ^3 E]- Below we will check directly that the terms containing 
curvatures in (j2)) are not zero on this solution. Hence, it is physically different from 
the solution of [TB] . 

For simplicity we concentrate on the r << r* region and do the following: we 
split h into two parts according to (fT^ . and for comparison with insist that tt 
has the form obtained in |16j. Using the results of |13j, for scales r <C r^., we get: 



where (3 = 3/2- 2{VS - 1) ~ 0.04, and a ~ ±0.84. While the vr field takes the 



Note that both hoo and tt given above are finite in the limit pOlll|) . It is the non- 
trivial part ()29|) that differentiates the solution of ^3 1^ from that of JHl- In 
particular, the curvatures R in Eq. Q are nonzero on this solution. 

The solution of [T8t IT^ has a number of interesting physical properties. Here 
we'd like to emphasize the effect of the gravitational mass screening [TH|ll4j. In the 
region r <^ r* the solution recovers the results of the Einstein theory with tiny, but 
potentially measurable deviations ^HIEIEI- Gravitational screening, on the other 
hand, manifests itself for r >> r*. In this region, the solution (Newton's potential) 
scales as ~ Tj.rM/f'^- Naively this seems counterintuitive, since it contradicts the 1/r 
scaling expected from perturbation theory in the region r^. <^ r <^ Tc [1 . However, 
as was explained in [THl I14j. there is no actual contradiction since the perturbative 
approach does not take into account the effect of gravitational mass screening of 
the source. Unlike in conventional GR, a static source in this theory produces a 
non-zero scalar curvature even outside of the source |l3i|. This curvature extends 
to scales ~ r*. Hence, the source is surrounded by a huge halo of scalar curvature. 
This halo screens the bare mass of the source. It is easy to estimate that the 
screening mass should be of the order of the bare mass itself [^: A deviation from 
the conventional GR metric at r -C r,, scales as mc^ruT (we ignore small /3 here.) 
This can give rise to the scalar curvature that scales as mc^/rMT'^^"^ . The curvature 
extends approximately to distances r ~ r^,. Because of this the integrated curvature 
scales as vrL(.^frMT* ~ tm, and the "effective mass" due to this curvature can be 
estimated to be ruM'^^ ~ M, i.e., of the order of the bare mass M itself! 

Given the above estimate, there could be either complete or almost complete 
"gravitational mass screening" for the problem at hand. Which one is realized in 
DGP? As we will argue below, the solution found in ^3 E] suggests that the 
screening is complete from the 4D point of view and is incomplete from the 5D 
perspective. Let us discuss this in more detail. From a 4D perspective (i.e., from 
the point of view of the induced metric on the brane) the solution behaves as a 




(29) 



form: 




(30) 
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spherically-symmetric distributions of mass/energy of radius r ~ r^,, with the bare 
mass M placed in the center, and the screening halo surrounding it. The solution 
[T3] shows that the potential at r > r,, has no "monopole moment" , i.e., it has no 1/r 
scaling. Instead, it has the "dipole moment" that scales as ~ r^rul'T^ ■, which fits the 
interpretation of a potential due to a 4D spherically-symmetric mass distribution of 
size ~ and zero net mass^. Hence, the 4D screening should be complete. The 
picture is slightly different from the 5D perspective. In the definition of the 5D ADM 
mass the off-diagonal component of the 5D metric of IT^ is also entering (this 
component vanishes on the brane and that is why it does not contribute to the 4D 
ADM mass). Because of this component, there is no complete screening of the 5D 
ADM mass, and the potential has the 5D "monopole" component which scales as 
~ r^TM/r'^ and yields the 5D ADM mass ~ M^rM/rcY^^ ||13j. The potential due to 
an incompletely screened 5D "monopole" smoothly matches onto the 4D potential 
due to a 4D "dipole" . 

6 Discussions and outlook 

It looks like the present model works in subtle ways. Yet, there are a number of 
issues that still need to be understood. We will outline some of them below. 

(1). DGP vs. massive gravity. The results of the present work show that in 
the decoupling limit massive gravity and the DGP model behave very differently. 
It would be instructive to understand this in terms of the Lagrangians of the two 
theories. It is useful to start by recalling some details of the decoupling limit in 4D 
massive gravity jTTj. The most convenient way is to use the Stiikelberg method and 
complete the FP mass term to a reparametrization invariant form J7] . This can be 
done order-by-order in powers of the fields. Let us call the 4D massive gravitational 
field f^^, the corresponding Stiickelberg vector field and its longitudinal part 0. 
Due to the mass term the field acquires a kinetic mixing term with the / field. 
The quadratic part of the action can be diagonalized by a conformal transformation 
ffiu = ffiu + Vfj.i^'t'- As a result, the (p field acquires its own kinetic term as well 
as coupling to the trace of the stress-tensor. The only nonlinear interactions that 
survive in the decoupling limit are those of (p, all the nonlinear terms containing 
/ vanish jTTj. The analog of the f^^, field in DGP is (7^,^, and the analog of 
is A''^, as it can be read off 0. Besides these fields there is an additional field 
N. In the linearized theory it enters the action Q as a Lagrange multiplier and 
enforces a constraint that is consistent with the linearized Bianchi identities. A 
linear combination of this field with the longitudinal component of A"^ is an analog 
of the longitudinal component of the field of massive gravity - these modes 
acquire their own kinetic terms through mixing with the tensor fields. However, the 
similarities end here. At the nonlinear level the A^ field ceases to be a Lagrange 

^We thank Gia Dvali for pointing out the dipole analogy to us. Notice, also that there is no 4D 
Birkhoff's theorem in this case. 
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multiplier, but enters © algebraically. Hence, it can be integrated out explicitly. 
The resulting action is a functional of g^y and only, however, it is very different 
from the action of massive gravity expressed in terms of f^^ and A^. The former 
contains non-local interactions between g^^j and A*"^. Moreover, if expanded on a flat 
background, it stays non-local and does not produce quadratic mixing terms between 
them. Under the circumstances, a way to proceed is to solve first equations for a 
nontrivial background, then expand the non-local action in perturbations about the 
background, and only then take the limit. This is, essentially, what we have done 
in the present work, but in easier terms of the equations of motion. 

(2) . Analyticity and the unitarity hound. The tree- level one-particle exchange 
amplitude in DGP is nonlocal from 4D point of view as it contains terms (□ + 
rricyj —^)~^ i22|j. The amplitude has a branch-cut because of the square root. As 
a result, there is a way do define a contour in the complex plane so that the pole 
in the amplitude ends up being on the second Riemann sheet (see, e.g., ^5). The 
usual 4D dispersion relation can be written for this amplitude, and it does respect 
4D analyticity. Having this established, one can consistently take uic — > 0, as it is 
done in the decoupling limit. How about the non-linear amplitudes? Generically, 
those contain terms that are singular in the rric —>■ limit (e.g., different powers of 
1 / rricy/ —O) These are the terms that make the perturbative expansion to break 
down at the scale r*. The results of the present work suggest that it is unlikely 
that the limiting theory below r* can be thought of a local 4D theory of massless 
helicity-2, helicity-1 and helicity-0 states. However, in the regions and/or in the 
backgrounds where the conventional perturbative expansion is valid (i.e, where the 
above singular terms can be ignored), it is reasonable to expect that 4D analyticity 
is respected. 

What is more certain, however, is that the Froissart bound cannot be assumed 
to hold even for the limiting theory. The reason being that in the decoupling limit 
there are long-range tensor interactions present, and the Froissart saturation can 
only be taking place if such interactions were absent. Clearly, the present theory of 
gravity has no mass-gap. 

(3) . (Important) miscellanea. The issue of the UV completion and quantum 
consistency of the model need more detailed studies (see, e.g., |^). In this respect, 
it would be interesting to pursue furtherer the issue of the string theory realization 
of brane induced gravity, perhaps, along the hues of Refs. [^01 EH |32] . 

Finally, in the present work we have not discussed the selfaccelerated universe 
jHSlEll- On issues of consistency of this approach will be reported in 
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